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Abstract 

In this paper the relationship between the problem of constructing the 
ground state energy for the quantum quartic oscillator and the problem of 
computing mean eigenvalue of large positively definite random hermitean 
matrices is established. This relationship enables one to present several 
more or less closed expressions for the oscillator energy. One of such 
expressions is given in the form of simple recurrence relations derived by 
means of the method of orthogonal polynomials which is one of the basic 
tools in the theory of random matrices. 
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1 Introduction 



In recent years there has been an interest in quantum mechanical quasi- exactly 
solvable problems since they seem related to many recent developments in the- 
oretical physics 0, ||. 

It was also shown that the Schrodinger equation for the energy eigenvalues 
for generic, that is non- quasi- exactly solvable potentials, may be approximated 
by an infinite sequence of potentials each of which is quasi- exactly solvable. More 
specifically, the solutions of the bound state equation for the quartic anharmonic 
oscillator are related to a sequence where the potential is a even polynomial of 
degree six. This fact enabled one to formulate a regular iterative method for 
constructing the ground state energy of quartic oscillator as an analitic function 
of coupling constant Q . 

In the present paper, it will be shown that the above mentioned sequence 
of approximations is naturally rewritten in terms of the large N limit for a 
random matrix^. The analysis presented here still is not a practical approach 
to the quartic anharmonic oscillator. However the techniques to analyze large 
random matrices had impressive developments in the past few years, so this 
approch seems promising. 

After rewriting the problem in terms of random matrix theory, it appears 
that the ground state energy of the quartic anharmonic oscillator is exactly 
provided by the first few terms in the topological expansion of a different "po- 
tential", in a zero-dimensional field theory. In other words, the exact sum of 
Feynman graphs of a one-dimensional model, which provide ground state en- 
ergy of the quartic anharmonic oscillator, are exactly given by the Feynman 
graphs of low genus, of a different zero-dimensional model. This feature is quite 
unexpected and is the second motivation of the present paper. 



2 Quasi-exactly solvable sextic oscillator model 

Consider a one-dimensional quantum sextic oscillator model with hamiltonian 

Hm = -^ + [b^ - aUM + 3)]x^ + 2abx^ + a^x^ (2.1) 
ox'' 

in which a > and b are real parameters and M is an arbitrarily chosen non- 



negative integer. It is known that model (2.1) is quasi-exactly solvable. This 
means that for any given M the Schrodinger equation for it admits algebraic 
solutions only for some limited part of the spectrum. The corresponding (alge- 

*This is a particular case of a more general result, according to which, any quasi-exactly 
solvable model is equivalent to a certain random matrix model M. 
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braically calculable) energy levels are given by the formula 



Em = {4.M + l)b + &aJ2^^^ (2-2) 



in which , . . . , S.M are real numbers satisfying the system of M algebraic equa- 
tions 

M 



with i — 1,...M. It can be shown that system (2.3) has M + 1 solutions 
describing energy levels with numbers 0, 2, ... , 2M. If the number of a level is 
2K then K of quantities S,i are positive and M — K oi them are negative. In 
particular, for the ground state all the quantities S,i should be negative. 



3 Quartic oscillator as a limiting case of quasi- 
exactly solvable sextic oscillators 



Let us now demonstrate that if M tends to infinity, the model (2.1) is reduced 
to the model with quartic anharmonicity. In order to demonstrate this let us 
assume that parameters a and b are functions of M. We choose this dependence 
in the form 

a^^gir\g,M), b^gif{g,M), (3.1) 

where f{g, M) is a function of two variables g and M implicitly determined 
from the equation 

/3(g, M) - .g-i/(.g, M) - 2(M + i). (3.2) 

In this case we will have 

6^ - a(4M+ 1) = 1, 2a5==(7, (3.3) 

iJjv^ = - — + (1 - 2a)a;2 + gx^ + a'^x^ (3.4) 

with a being the root of the cubic equation 

^2 ^4a2[l + a(l + 4M)] (3.5) 
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in the interval < a < g/2. Since 



1 g2/3 

a fa — for Af —> oo, g fixed (3-6) 



2 (2M)i/3 



in the hmit M — > oo the model (2.1) reduces to the quartic oscillator model 
with hamiltonian 

H= lim HM^-— + x^+gx\ (3.7) 



Now what about the solutions of the model ( |3.7| )? It is quit e ob vio us th at 
at least all its even solutions can be reco vered from formulas ( |2.2| ) - (O) if 
we substitute in them the expressions (3.1) for a and b and then take the limit 
Af — > cxD. For the sake of simplicity we can restrict ourselves by discussing 
only the ground state solution for which all the numbers ^i, . . . should be 
negative. But before going over to the limit Af ^ oo, it is desirable to make 
(for the sake of further convenience) some appropriate rescaling of quantities 
^1) • ■ ■ I S,M^ introducing instead of them the new quantities 771, ... , rjM by the 
formulas 

C« =-5"*/'(3,Af)r7„ z = l,...,Ar. (3.8) 



After this the expression for the ground state energy in the model (3.7) takes 
the form 



E= \\m Em = {AM + \)g^f{g,M)-[l- SM{ni.---.nM)]. (3.9) 

M— >oo Af^oo 



where 



M 



i=l I 

and 771 , ... , f]M are positive solutions of the system 

M 

E + — + /3(<?,A/)(l-77,)=0, (3.11) 

with 1 = 1,.. .M. 

Taking into account the facts that the large M asymptotics of function 
f[9,M) is 

/(5,Af) = 2i(Af + 1/4)3 (1 + \ ^ + ...1, (3.12) 

^ ' ^ ' ' \ 3[2g(M + l/4)]i J ^ ^ 
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and the ground state energy E of the quartic anharmonic osciUator is a finite 
number, we can conclude that for large values of M the quantity S'm(^i, ■ • ■ , Vm) 
should behave as 

E 1 

Sm{vu . . . ,?7m) = 1 - —-—T ■ —pr + ■■■ (3-13) 

4(25)3 Ms 

In other words, if 5^/(^1, ■•■ , Vm) is already known and we would like to find the 
ground state energy of the quartic anharmonic oscillator, it is sufficient to find 
the first non- vanishing correction to unity in the expansion of Sm ivii ■ ■ ■ ,Vm) 
in negative powers of M. All higher corrections will be irrelevant to us. 



Integral representation for the ground state 
energy in quartic oscillator model 



The resu lt o f the previous section enables one to replace S'm(??i, • ■ • , ?7a/) 
formula (|3.9| ) by another quantity which we denote by Sm and define as 



M 

E 



<Vi> _ JT Dt]Sm{vi, Vm) exp[-V{T]i, tjm)] 



(4.1) 



where 



and 



^M + 1/4 j;^ D7^exp[~V{m,---,m)] 

M 

Df]^Y[ dTj, (4.2) 



i=l 



M , A/ M 



1^(771, . . . , t^m) = - In - %| - ^ X] 1^ l^'l + M) f^im - 1)' (4.3) 

i^k i—1 i—1 

The reason for such a replacement can be explaned as follows. It is clear that 



the exponent (4.3) has a minimum which is saturated by numbers rj satisfying 



the system of equations (3.11). As follows from these equations, all the numbers 
fj sh ould be of order of unity, and therefore, the most relevant values of exponent 
(|J) should be of order M^. This means that for large values of M the multiple 
integrals in (^^) can be evaluated by the saddle point method and the correction 
to the zero order approximation should be of order M~^. Now note that the 



zero order approximation for (4.1) is nothing else than 5j\/(^i, . . . ,?7m)- The 
fact that the difference between SmIvi, ■ ■ ■ ,Vm) and Sm is of order (i.e. 
less than M~'s ) makes it possible to replace Sm{vit--tVm) by Sm in the 
expression for E. This enables one to write the following simple formula for the 
ground state energy in the quartic oscillator model (|3.7|): 



E = lim (4M + l)g3 /(g, M) • [1 - Sm] , (4.4) 

M— »co 



where Sm is defined by formulas (4.1) - (4.3) 
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5 Random matrix model 



In this section we recall that the expression (4.1) for Sm naturally arises in the 
theory of random matrix models. 

Consider a zero-dimensional "field-theoretical" model in which the role of a 
field is played by the random M x M positively definite hermitean matrix A. 
Let us choose the action of the model in the form 

WM[A]^TT[^^\nA + fig,M)iA-lf] (5.1) 

which is obviously invariant under unitary transformations A — *■ UAU"^. The 
simplest invariant quantities in this model are the eigenvalues of the matrix A. 
Consider the problem of computing A ~ (Tr A) / (M + 1/4) . The result can be 
represented in the form of a multiple integral 

j_ 1 J DATTAeM-WM[A]} 

M-fl/4 f DAexp{-WM[A]} ^'^ 

which can be simplified if one uses the well-known Faddeev-Popov method for 
extracting from (|]^) the volume of the unitary group. Using formulas 

A^U ■diag{m,---,mi}-U-' (5.3) 

and 

M 

DA = Y[dT],l[{rj,-7jkfDU (5.4) 

i—l i<k 

in which r]i are the eigenvalues of A and DU denotes the measure on the unitary 
group, we can easily reduce the expression ( p3.2[ ) to the form 

^ _ /o°° DvSm{vi, Vm) exp[-V{i]i, tjm)] 
Drjexpl-Vif]!, . . . ,?7m)] 



i n wh ich D rj, S m{vii ■ ■ ■ , Vm) and V^iji, . . . , t/m) are given by formulas (4.2), 
(3.1C) and (4.3), respectively. We see that formulas ( |5.5D and ( |4.1| ) exactly 



coincide! This means that the problem of finding the ground state energy in the 



quartic oscillator model (3.7) is equivalent to the problem of finding the mean 



eigenvalue of the random matrix A in the theory with the action (5.1) 



6 Orthogonal polynomials and recurrence rela- 
tions 



The equivalence of models (3.7) and ( ^.l| ) enables one to simplify the expression 
(4.4) by using the method of orthogonal polynomials which is one of the basic 
tools in the theory of random matrix models [Q. 
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The idea of this method can be formulated as follows. Let us rewrite the 



expression (4.1) in the form 

g ^ /o°° UfLi dri.SMjrji, ijm) Ilf<kiv^ ~ Vkf exp[-K)(r/i, . . . , ??m)] 
/o^ n,=i drl^ Y{T<ki'n^ - VkY exp[-yo('7i, ■ • ■ , '7m)] 

with 

M M 

i 

2" 



M M 

Voivu . . . , ??Af) = ^ ^ In h\ + f{g, M) J2{V, - 1)' (6.2) 



i=l 

and note that 



l[iv^-Vky = [det||7?f-i||,,fc=i,...,M]' = [det||Pfc_i(77,)||,,,.=i,...,A/]' (6.3) 



i<k 

where Pniv) denote arbitrary linearly independent polynomials of degrees n 
normalized as P„(77) = + . . .. Because of the arbitrariness of these polynomials 
we can consider them (without any loss of generality) as polynomials orthogonal 
with the weight |?7|^/^ exp[— /^((7, M)(77 — 1)^], i.e. requiring the condition 

i 

dvP^{rj)P^{v)\rj\^/^exp[~f{g,M){rj-l)^] = h^Smn (6.4) 

This immediately leads us to the formula 

. 1 J^°°dr,^Pl{^)\r,\^/^eM-f{g,M)i^-m 

= wnr^ 2: (6.5) 

' m— 

By using well known methods, [H, recurrence relations are derived for the 



orthogonal polynomials, which allow the evaluation of the expression (3.£). The 
polynomials are given by a three terms equation 



rjPniv) - Pn+liv) + Mv)Pniv) + Rn{v)Pn-liv) (6-6) 

with Rn{r]) — j^^- It follows that 

M-l 

' k=n 

It is also customary to introduce the states \n > 

< ri\n >= -P„(?7), where < m\n >= S^n (6-8) 
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One obtains a local string equation 

< n\r]W'{r])\n >= 2n + 1 

that is 

2 fig, M)[Al -A„ + Rn+i + i?„] = 2n + ^ 
and a non-local equation 



(6.9) 
(6.10) 



J^^k = VRn+i <n + l\rjW'i7j)\n 



fe=0 



- i?„+i2/3(5,M)(A„+i+A.-l) 
which leads to a further expression for Sm 



-Rm[Am + Am-1 - 1] 



(6.11) 



(6.12) 



Af + 1/4 

If the coefficients Ak, Rk are known up to fc = ri, eq. ( 3.10| ) allows the 
determination of Rn+i , next A„+i may be evaluated by eq. ( 6.11 ). The lowest 
values are 



An 



Ml 

Mo 



Ai = l 



where jin are the moments 

M„(5,M)= / dxx"+i/2e-/'(3^^'^)(^-i) 



i?0 = 



Mo 



(6.13) 



(6.14) 



Since 5'a/ in eq.(6.12) determines the ground energy E{g) of the quartic an- 
harmonic oscillator (4.4) for large values of Af, the above iterations do not 
seem very useful. It might seem more promising to derive a recurrence relation 
immediately for the quantities 



One easily obtains: 

{qn+l - qn ) [qn+l - qn ~ 1) + 



(6.15) 



1 



q-ii+i 



9n+l " 9n-l — 1 



2p{g,Myq^+2-qn-l 
- 2n- 3/2] = 



(6.16) 
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where qn+2 is determined in terms of qn+i, Qn, Qn-i- The lowest values are 
qo = 0, = Ao, q2 = Ao + Ai 

Thus the final formula for the ground state energy of the quartic oscillator 
model (3.7) has the form 

E= lim Em = lim 4(2M''g)* 



M- 



M- 



1 - 



qM 



M+l/A 

where qM can be found from the recurrence relation ( 6.16| ). 



(6.17) 



7 Conclusion 



The equations ( |6.16 ) and ( ^.17 ) obtained in the above Section are exact and, 
in principle, provide an exact evaluation of E(g) (i.e. with an arbitrarily high 
accuracy). It seems, however, that these equations have yet a purely theoreti- 
cal significance, because the practical use of these formulas encounters serious 
difficulties. 

The first difficulty follows from the fact that the deviation of the exact 
value of E from its approximations Em obtained for finite M's is of order 
Indeed, remember, that the energy E appeared as a coefficient in front 
of M""'/^ in the expansion of Sm{vi, ■ ■ ■ ,fjM) (see formula (3.13) ), and we 
replaced Sm{vi> ■ ■ ■ tVm) by Sm after neglecting corrections of order (see 
Section 4). This means that the sequence of Ei\i converges to E very slowly. 
For example, if we want to compute E with 1 percent accuracy, we should take 
M = 1000, and this means that we need 1000 steps of the recurrence procedure 
( |6l^ ). The recurrence relations (|6.16[) are, however, rather unstable, and the 



iteration of q„ from low to high values of n does not seem practical, because of 
the rapid accumulation of errors. 

The second difficulty lies in the fact that the direct analysis of the asymptotic 
behaviour of the difference equation ( 6.16 ) is non trivial. It is rather easy to 
obtain the leading behaviour of numbers g„ for large n, considering g„ as an 
analytic function of n, qn — q{n) and solving the associated differential equation 
for q{n). The result, g„ ~ n-|-l/4, does not contradict formula ( 6.17 ). However, 
a naive attempt of findi ng in such a way the non-leading behaviour of q{n)^ which 
is relevant for eq.( 6.17 ), leads to a priori false result. 

It seems that both above mentioned difficulties are closely related to each 
other. It is possible that function q{n) is continuous, but not differentiable. 
In other words, the first (or, perhaps, the second) derivative of q{n) might 
be a random function of n. In this case, the approximation of the recurrence 



relations (3.1(;) by a differential equation is, generally, rather rough, and cannot 
be used for studying any subtle properties of function q„. If the behaviour of the 
differences q„ — (the derivatives of q{n)) is really stochastic, then this fact 
would easily explain the rapid accumulation of errors in numerical computations 
of the ground state energy of quartic oscillator by formulas (^.16|) and (3.17). 
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In conclusion of this Section we demonstrate that the problem of determi- 
nation of a non-leading behaviour of numbers qn is related to the problem of 
finding the form of the second local string equation for the recurrence relations. 
In order to show this, let us first provide an approximate evaluation of Sm 
for large M by means of the saddle point method. At leading order, for large 
M, one may replace the discrete eigenvalues rji with the continuous function 
r]{i/M) = rji and the algebraic system ( 3.11) with the singular integral equation 



^2/Z7:::^^7rV + TTTTT^ + ^^7^ [1 - ^(^)] = (7-1) 



£ T]{x) - r]{y) ' 4M7j{x) ' M 

drj 



which, after the introduction of the eigenvalue density p{r}) — becomes 



This is the saddle-point approximation for the large-M limit of the matrix model 
(pn|). Its solution with < A < B was described in together with the 
evaluation of expectation values. One obtains, for large AI 

J2<m>-M r ^,p{^,)d^Ji = {B~Af{A + B-l)^-^^^^ (7.3) 
i=l •'^ 

where the extrema of the support are solution of the algebraic equations 

f3 



2VABiA + B -2)f-\g,M) = 1 (7.4) 



From the system (7.4) one easily evaluates the large M expansions 



A 



Af2/3 Af4/3 



^ = 2 + ]i^ + ]i^ + - (7-5) 

To estimate the present saddle-point approximation, we quote the first coeffi- 
cients in the limit of large values of g 

-2 = y^^[l + 0{g-^^')] , ^'2 = ^^[l + 0(5-^/^)] (7.6) 



Then 



M 1 Q 1 

2 — 1 ' 

E{g) ^ gy^h^l + 0{g-^'^)] (7.7) 
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The comparison of the approximate result (7.7) with the known first term of 
the strong couphng expansion ||] 

E{g) = 22/3gi/3[o.667+ 0(g-2/3)] (7,8) 

shows that the fractional powers of the coupling constant g are correctly repro- 
duced, but the relative error for the first coefficient in the strong coupling limit 
of the ground state energy is about 32%. This fact is not surprizing because we 
started with a continuous version of the saddle point equation, which is only an 
approximation to the true discrete one. 

Let us now consider the equivalent approximation by recurrence relation. If 
we neglect the contributions in integration by parts at ?7 = 0, discussed in a 
second local string equation is obtained Q : 

<n\W'{ri)\n>=Q (7.9) 

and evaluating it in the continuum approximation ( that is A„ = yl(-||-) , i?„ = 
R{^)) and neglecting the derivatives, it reads 

4/3(5, M)[A{x) - 1] - I (7.10) 

The same approximation on the first local string eq.( 6.10| ) is 

'^-l^^[A'{x) - A{x) + 2R{x)] - 2x + A (7.II) 

This algebraic system, for large values n ^ M that is a; ~ 1 reproduces the 
system (7.4) after the usual correspondence 



16 

^(1) ^ ^ (7. 12) 

From the above exercise we learn that to improve on the approximate evaluation 
( [7.7] ) by use of the recurrence relations, it is likely to need an improved evaluation 
of the second string equation, by including the derivatives and the boundary 
contributions at ?/ — 0. We intend to consider this question in fortcoming 
publications. 
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